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ABSTRACT 


This  computer  program  provides  launch  to  burnout,  two- 
dimensional,  point  mass  trajectories.  It  is  a  very  flexible  program 
in  that  it  may  be  based  on  a  large  variety  of  conditions  used  to  control 
the  program.  The  standard  atmospheric  conditions  and  gravitational 
field  are  included. 

The  two  versions  of  this  program  are  written  in  FORTRAN  II. 
The  most  flexible  version.  SP11A,  is  card  tape  approach,  but  it 
requires  a  longer  "run"  time.  The  SP11B  is  the  all-card  approach. 

It  can  be  used  on  either  the  IBM  1410  computer  or  the  IBM  1620 
equipped  to  use  FORTRAN  II. 
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INTRODUCTION 


This  computer  program  was  developed  as  a  tool  for  studying  boost 
phase  trajectories.  It  is  used  to  investigate  the  interaction  of  the  basic 
ballistic  parameters  such  as  thrust,  mass,  drag,  gravity,  and  wind. 

The  present  program  allows  only  two  degrees  of  freedom  of  a  point 
mass.  However,  a  subsequent  program  will  provide  for  four  degrees 
of  freedom. 

The  output  of  this  program  includes  time  of  flight,  range, 
altitude,  acceleration,  velocity,  flight  path  angle  and  other  missile 
and  trajectory-related  variables.  Because  of  the  relatively  short 
distance  traveled  and  the  short  flight  times  associated  with  the  launch 
phase,  it  was  not  considered  necessary  to  use  an  earth-fixed  geocentric 
co-ordinate  system.  The  two-dimensional  co-ordinate  system  used 
in  the  program  has  its  origin  at  the  launch  point.  The  X-axis  represents 
range  and  is  tangent  to  the  earth's  surface  at  the  launch  site.  The  Y-axis 
represents  altitude  and  is  normal  to  the  earth's  surface  at  the  launch 
site. 
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ANALYSIS 


Initial  Conditions 

Almost  any  combination  of  initial  conditions  can  be  used.  In 
addition  to  standard  trajectories  based  on  surface -launched  missiles 
having  zero  initial  velocities,  this  program  can  be  used  to  calculate 
second  stage  boost  trajectories  or  air-launched  missiles  and  other 
trajectories  based  on  non-standard  conditions. 

The  items  of  information  to  be  punched  on  IBM  cards  are: 


Items 

RMAS 

BIWT 

FWT 

VL 

AL 

A 

TIME 

DELT 

D 

THETAO 


Comments 

Propellant  weight  flow  rate  before  the  thrust 
reduction  period  -  lb/sec 

Initial  mass  -  slugs 

Missile  weight  at  burnout  -  lb 

Initial  longitudinal  velocity  -  ft/sec 

Missile  length  -  ft 

Vehicle  Reference  Area  -  ft* 

Time  thrust  reduction  begins  -  sec 

Computation  interval  -  sec 

Diameter  of  missile  •  ft 

Angle  of  the  velocity  vector  with  respect  to 
the  horizontal  measured  counter-clockwise 
at  launch  -  radians 
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Items 


Comments 


AK1 

AK4 

THRUST 
ALTO 
DRAGO 
A  LAM 
G 

ALAMI 

TP 

FTIME 


ALA 

FRMAS 

SSI 


SS2 


Angle  of  the  velocity  vector  with  respect  to  the 
horizontal  measured  counter-clockwise  at 
burnout  -  radians 

Pitch  angle  to  be  imposed  on  the  missile  after 
time  (TP)  -  radians 

Initial  thrust  -  lb 

Launch  altitude  above  earth's  surface  -  ft 
Initial  drag  -  lb 

First  limit  on  drag  coefficient  curve 

Acceleration  due  to  gravity  at  launch  altitude  - 
ft /sec2 

Second  limit  on  drag  coefficient  curve 

Time  to  impose  the  pitch  angle  AK4 

Time  of  burnout  -  sec  (If  no  thrust  reduction 
period  is  desired,  let  FTIME  and  TIME  be 
the  Barne. ) 

The  longitudinal  acceleration  to  be  imposed 
on  the  missile  during  the  thrust  reduction 
period  -  ft/ sec2 

Propellant  weight  flow  rate  during  the  thrust 
reduction  period  -  lb/sec 

Control  switch.  If  SSI  >  0,  use  a  thrust 
curve  as  f(lime)  on  tape  drive  6.  If  SS1>0, 
either  a  constant  thrust  or  a  computed  thrust 
will  be  used.  For  this  decision  check  THRU 
and  TTTT. 

Control  switch.  If  SS2  >  0,  use  a  wind  profile 
curve  as  a  f(aititude)  on  tape  drive  5.  If  SS2  <  0, 
no  wind  consideration  will  be  made. 
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Items 


Comments 


553  Control  switch.  If  SS3  >  0,  use  a  pitch  program 

with  changes  in  6  as  a  f(time)  on  tape  drive  7. 

If  SS3  <0,0  will  be  computed. 

554  Control  switch.  If  SS4  >  0,  use  an  acceleration 

curve  as  a  f(time)  on  tape  drive  8.  If  SS4  <  0, 
the  acceleration  will  be  calculated. 

555  Control  switch.  If  SS5  >  0,  time,  A^,.  Vj^, 

Y  and  X  will  be  punched  in  a  card  for  each 
computation  interval. 


For  the  version  of  the  program  which  utilizes  tapes  as  well  as 
cards,  the  items  of  information  to  be  stored  on  magnetic  tapes  are  as 
follows: 

1.  Drag  coefficient  as  a  function  of  Mach  number 

2.  Wind  velocity  profile  as  a  function  of  altitude 

3.  Acceleration  as  a  function  of  time 

4.  Thrust  as  a  function  of  time 

5.  Flight  path  angle  as  a  function  of  time 
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Development  of  Equations 


1.  Thruet 

(a)  When  thrust  is  used  to  drive  the  program,  it  is  used  either  as 
a  constant  or  a  tabulated  variable.  Input  into  the  program  is 


based  on 

the 

following  relationship: 

A-p 

= 

T/m„ 

where: 

Aip 

s 

Acceleration  due  to  thrust 

T 

= 

Thrust 

mn 

s 

Mass  at  any  time 

(b)  When  a  pre-programmed  longitudinal  acceleration  is  used  to 
drive  the  program,  thrust  becomes  an  output  and  is  described 
by  equation  (2): 

e 

T  *  ^2 — Ul!  (Ai  +  An  +  g  sin  0  -  Au»  cos  0)  (2) 

cos  o 

where: 

m0  «  Initial  mass 
m  ■  Mass  flow  rate 
t  *  Burn  time 
o  *  Thrust  alignment 

■  Longitudinal  acceleration 
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Ad  ■  Deceleration  due  to  drag 

g  *  Acceleration  due  to  local  gravity 

6  ■  Flight  path  angle 

Ayf  =  Acceleration  due  to  wind 


1 

1 

1 

1 

I 

I 

I 
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2.  Drag 

Drag  ie  introduced  into  the  program  by  means  of  the  deceleration 
equation  below: 

AD  *  ipV^CpA  (J) 

mn 

where: 

p  s  local  air  density 

*  longitudinal  velocity 
Cq  -  drag  coefficient  (f{Ma) 

A  =  reference  area 

The  drag  coefficient  function  is  based  on  a  sero  angle  of  attack. 

It  is  a  function  of  Mach  number  and  is  included  in  the  program 
on  tape  or  as  a  subroutine  depending  on  the  program  version. 

3.  Flight  Path 

The  program  is  set  up  in  such  a  way  that  flight  path  can  be  an 
input  as  a  function  of  altitude  or  an  output  based  on  equation  (4): 

6  =  tan’1  Vy/Vx  (4) 
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4. 


Wind  Effects 


Accelerations  which  are  a  function  of  wind  are  included  in  the 
program  by  way  of  the  relationship  described  below: 


_ 


m. 


where: 

P  Vw2 

Pj-j  =  dynamic  pressure  =  — —— 
Ag  =  exposed  area  =  LDsin  0 
Vyy  =  velocity  of  the  wind 
L  =  length  of  missile 

D  =  diameter  of  the  missile 


(5) 


5.  Composite  Accelerations 

Summations  of  the  accelerations  along  the  X  and  Y  co-ordinates 
are  obtained  from  equations  (6)  and  (7). 

A\  -  A>j>  cos  (0  -  o)  -  Aq  cos  0  +  (6) 

Ay  s  A-p  sin  (0  -  a)  -  Aj^  sin  0  -  g  (7) 

6.  Longitudinal  Acceleration 

Longitudinal  Acceleration  is  then  a  resultant  of  the  X  and  Y 
accelerations  as  described  by  equation  (8): 


Al  =  Ax  cos  0  +  Ay  sin  0 


(8) 
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7.  Velocities 


Velocities  are 

accelerations. 

obtained  from  integration  of  the  various 

t 

vx  « 

/  AX  dt 
o 

t 

(9) 

ii 

> 

> 

/  ay  dt 

O 

t 

(10) 

VL  = 

i  al  dl 

(ID 

o 


8.  Distance  Traveled 

The  X  and  Y  co-ordinates  and  the  longitudinal  distance  traveled 
are  obtained  by  integration  of  the  corresponding  velocities. 


t 

X  = 

/vxd. 

o 

t 

(12) 

Y  * 

/  vv-“ 
o 

t 

(13) 

L  = 

f  Vudt 
o 

(M) 
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FORTRAN  STATEMENTS 
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CARO  TO  PRINT  80/80 
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FORTRAN  RUN 

HOP  16 

LAUNCH  PHASE  TRAJECTORY  BY  W  B  WARREN  FOR  RESEARCH  STAFF 

24  AUGUST  1962 


PROGRAM  NO.  649 

RMAS  »  PROPELLANT  MASS  RATE  OF  FLOW  -  LBS/SEC 
BIWT  H  TOTAL  INITIAL  WEIGHT  -POUNDS 
FWT  »  MISSILE  WEIGHT  AT  BURNOUT  -POUNDS 
VL  0  INITIAL  LONGITUDINAL  VELOCITY 

AL  0  MISSILE  LENGTH  -FEET 

A  U  VEHICLE  REFERENCE  AREA  -FEET 

TIME  0  TIME  OF  BURNOUT  -  SECONDS 

C  0  OIAMETER  OF  MISSILE  -  FEET 

OTHETA  0  ANGLE  OF  THE  VELOCITY  VECTOR  WITH  RESPECT  TO  THE  HORI/ONTA 
1L  -RAOIANS 

AK1  0  THETA  AT  BURNOUT  -  RADIANS 
ALPHA  0  ANGLE  OF  ATTACK 
THRUST  0  POUNDS  FORCE 
ALTC  0  LAUNCH  ALTITUOE 
ORAGC  0  INITIAL  DRAG 

ALAM  0  LIMIT  ON  DRAG  COEFFICIENT  CURVE 
ALAMl  0  LIMIT  ON  DRAG  COEFFICIENT  CURVE 
G  0  GRAVITATIONAL  FORCE  FOR  LAUNCH  ALTITUDE 
OUNITS  OF  TABLES  USEC  -SPELD  OF  SOlJND-F  T/SEC  ,  DENSI TY-LB/F T**3 ,  GRA 
IVITY  -FT/SEC*«2,  ALTITUOE  -FI ET 


REAL  100, RMAS, BIWT, FWT, VL , AL 

RE  AC  100 ,A, TIME, 0ELT,D,THCTA 

RE AC  lOO , AK 1 , AK4, THRUST , AL TO, DR AGO 

REAL  1 00, AL AM, G, ALAM l, TP, FT  IME 

RE  AC  100 , AL A, FRMAS 

REAL  100,SS1,SS2,SS3,SS4,SS5 

PRINT  1000, RMAS, BIWT, FWT, VL.AL 

PRIM  1001,A,TIME,0£LT,D,  fHLTA 

PRINT  1007, AK1, AK 4, THRUST, AL  TO, DR AGO 

1000  FORMAT  X 1H1 ,6HRM AS  , E 1 5 . 8, ?X , 6H I WT  , E l 5. 8, 2X ,6HF WT  »E15.8»2X, 

16HVL  ,E15.R,2X,6H  L  .ElS.Hn 

1001  FORMAT  «1H  ,6H  A  , E l 5. B, ?X , 6HT IME  ,E 15. 8*2X,6H0ELT  ,EIS.8,2X, 

16HDIAH  ,E15.8, 2X.6HTHETA  ,EIS.8t) 

1002  FORMAT  tlH  ,«H  Kl  ,E15.8,2X,6M  K4  ,E l 5. 8, 2X,6HTHRUST ,E IS. R , 2 X , 

16HALT0  ,E15.8,2X,6HDKAGO  ,Ll5.0a 

SENSE  LIGHT  2 
XXXX  IT  0. 

THCTO  ITTHETA 
AK'J  0  30. 

ALPHA  0  0. 

100  TCRMAI  *5E15.8a 
T  i  U. 

CCCC  f  0. 

BILL  0  0. 
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CARO  TO  PRINT  80/80 


IPAGE  A  0 
CCNTR  A  0. 

SENSE  LIGHT  4 
AP  ASS  A  BIWT 
ALT  A  ALTO 

1400  RE AC  INPUT  TAPE  l, 2000, AL T 1 , DEN  1, GRAV 1 , THRUl , SOSl 

2000  FCRPAT  S5El5.8n 

I FSALT 1-500000 .048 1,481,482 
482  SCSIAO. 

481  DEN1ACEN1*. 0310807453 
IF  SALT  —ALT  la  1 , 2, 3 

3  ALT2AALT1 
0EN2 ADEN  1 
GRAV2AGR AVI 
THRU2  A  THRUl 
SCS2  A  S0S1 
GC  TC  1400 

2  ALT  A  ALT  1 
GRAV  AGRAV1 
DEN  ACEN1 
ALT2  A  ALT  1 
GRAV2  A  GRAV 1 
CEN2  A  CEN1 
SCS2  A  SCSI 
SCS  ASCS1 
40  IF  XSSia7,7«6 

l  ALRAT  A  *ALT-ALT2o/tALT 1-ALT2D 
OEN  A  0EN2  -  * ALRAT *XCEN2-0£N Iob 
GRAV A  GRAV2-XALRAT*XGRAV2-GRAVlcu 
IF  SSCS1  -  SOS 20  1301,1302,  1303 

1302  SCS  A  SCSI 
GC  TC  1304 

1301  SCS  A  SCS2  -  X ALRAT  •  XS0S2  -  SUSIcd 
GC  TC  1304 

1303  SCS  A  SCS2  €  X ALRAT  •  XSOSl  -  S0S2nc 

1304  BACKSPACE  1 
GC  TC  40 

7  IF  STHRUSTo  8,9,10 

9  SENSE  LIGHT  l 

GC  TC  23 

8  GC  TC  10 

10  THRUPTHRUST 
GC  1C  23 

6  RE AC  INPUT  TAPE  6, 2 100, T l , THRUl 

2100  FCRPAT  X2E15.Ro 

IF  X1-TIO7200, 7201,7202 

7201  THRU  A  THRUl 
THRU?  A  THRUl 
GC  TC  23 

7200  THRU2ATEKU1 
T2  A  T 1 
GC  1C  6 

7202  TRAT  AXT-T2U/XT l-T2o 
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CARD  ro  PRINT  80/80 


THRU  H  THRU 2  £  XXTHRU2-THRU In* TRATn 
BACKSPACE  6 
GC  TC  23 
23  AMVL/SCS 

lF*AMn87,76, 11 
76  ORAGH/O. 

GC  TC  151 

11  REAC  INPUT  TAPE  2,  2200, AM  1, DRAG  1 
IF  UP  -  Al AM. 15,  15,14 

14  IF  UM- AL  AM  lnl6, 16, 15 

15  IF  UM-AMinl7,  18,  19 

19  0RAG2  »  DRAG  1 
AM2  *  AMI 

GC  TC  11 
18  DRAG  *  CRAG l 
AMACM  H  AMI 
BACKSPACE  2 
GC  TC  151 

17  AMRAT  0  T!AM-AM2o/t AMI-AM2U 

DRAG  #  CRAG 2- AMRAT *t0RAG2-0K AG  1 o 
AMACH  0  AM 
BACKSPACE  2 
GC  TC  151 

16  IF  7  AM-  AM1o20,21, 22 

20  AMRAT  *  '{ AM- AM?o/XAM  1— AM2n 

CK AC  0  CRAG2£tAMRAT«SCRAGl-URAG?DD 
AMACH  0  AM 
BACKSPACE  2 
GC  1C  151 

21  UH AC  0  CRAG  1 
AMACH  0  AMI 
GC  TC  151 

22  CP AC 2  A  0KAG1 
AM?  «  AMI 

GC  1C.  11 

151  IF  «.SS?n?7 ,27*26 

27  CPW  *0. 

GC  TC  2222 

26  RF Al  INPUT  TAPE  5 , 2200, HE IG I , VW l 
IF  <AlT-HUGla?8,?9,30 

30  HE  I G?  0  MTIG1 
VW?  H  VW  1 

GC  TC  26 

28  HRAT  0  ‘SHClGl-ALTa/tMEIGl-HE  lG2c 
VW  0  VW l-«MRATa*XVW l-VW2o 

31  CPW  mCENO*VW**2nn/2. 

BACKSPACE  5 

GC  TC  222? 

29  VW  0  VWl 
GC  H.  31 

2222  AC#VU./2.o«CEM*XVL**2B«ORAG*Ao/AMASS 
IF  IT  -  TPal53, 153,915 
915  IF  INS  30/77,7  77, 32 
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CARD  TO  PRINT  80/80 


777  IF  XALPHAU152,  153,  153 
32  RE  AC  INPUT  TAPE  7, 2300# T  IM1 , THE T l 
2300  FCRMT  X2E15.8B 

IFXT-TIPt£>34,35,36 

35  THEI2  A  THET 1 
GC  TC  152 

36  TIH2  # T I M I 
THET2ATHET l 
GC  TC  32 

34  TUATAXT-TIM2a/*TlMl-r  lM2n 

THET2  ATHET2-XTRAT*XTHET2-THET  UB 

152  RE AC  INPUT  TAPE  9, 2 300, T IM  1, ALPH 1 
IF  *T-TIM3o37,38,39 

38  ALPFA  AALPH1 
GO  TC  153 

39  T I P4AT IP3 
ALPF2  A ALPH1 
GC  TC  152 

37  T3RAT  AXT-T  1M40/XT  IH  J-T  IM4a 

ALPHA  U  ALPH2  -XT 3R A T *X AL PH2- AL PH Iud 

153  AE  A  AL*C»S INFXTHET An 
AWAXCPW* AEn/ APASS 

IF  XSENSE  LIGHT  lo  700,701 
700  IF  *SS4a  702,702,703 

703  Rt AC  INPUT  TAPE  8 , 2200, F T IME l , AL A  1 
2200  FGRHAT  X2E15.80 

IF  * T-F T IMEO704.705,  706 
706  FTIHE2  AFT  I  ME  1 
ALA2  A  ALA l 
GO  IC  703 
705  ALA  AALA1 
AL A2A ALA  1 
FT IHE2AFT IME I 
T  AFTIMEl 
GC  TC  702 

704  TIMER  A  XT-FT IME2o/XFTIMEl  -  FT IME  2n 
ALA  »ALA2  -XT IMER*XFT IME2  -FIlMElon 
BACKSPACE  8 

702  IF  XCCCCo769,769,770 

769  REAC  100, THRU, TTTT 
CCCC  A  CCCC  t  1. 

GC  TC  911 

770  IF  )TTTTa501,501,91l 

501  THRUA  XAMASS/COSFXALPHAoo*XAL  AEADCAG*  SINFXTHE  TAa-AW*COSFXTHE  TAau 
911  AT AlHRU/AMASS 

AX  A  AT •CUSFXTHET A  -  ALPHAa  -  XAD»COSFXTHE  TAanC  AW 
IF  X AX  -  .10925,925,926 

925  AX  A  0. 

926  AY  A  AT *S INFXTHET  A  -  ALPHAo  -  XAO*SINF  (THE  TAaa-GRA V 
BL  A  AX*  COSFXTHETAa  £  AY*  SINFXTHETAo 

IF  XSENSE  LIGHT  4051,52 
51  VX  AO. 

VY  A  0. 
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B  V  L  M  BL 
X  *  0. 

Y  *  0. 

CVL  «  0. 

SENSE  LIGHT  4 
GC  TC  54 

52  VX  #VX  C  SAX  •  CELT  □ 

VY  0  VY  £  SAY  •  CELT n 
902  BYL  H  SCRTFSVX*«2  £  VY»*2o 
X  0  X  £  VX*CELT 

Y  *  Y  £  VY  •  CELT 

CVL  X  SCRTFSX«»2  £  Y*«2o 
54  AG  0  GRAV 

296  CCKTR  ACCNTR  £  I. 

IF  *CCNTR-l.o297, 298,299 

297  PAUSE 

299  IF  SCCNTR  -7.0399, 399, 365 
365  CCNIR  0  0. 

GC  IC  296 

298  1  PACE  0  I  PAGE  £  1 
PRIM  5000,  IPAGE 

5000  FORMT  SlHl ,  1 10X,  5  HP  AGE  ,  I  3o 
399  PRINT  4999 
4999  FCRPAT  SlH  ,//0 
5701  A A A A I  0  ALPHA/ .01745 
AAA A?  0  THETA/. 01745 
PRIM  5001, T,X,Y, CVL,  AAAA1 

50010FGRMT  SlH  ,6HT  IME  ,E15.8,?X,3H  X  ,EIS.B,?X,3H  Y  ,E15.8,2X,3H  L  , 
IE  15.8, 2X,6H ALPHA  ,EI5.8a 
ANN  0  X/6080. 

ALT  0  Y  £  ALTO 

PRIM  5002, THRU, VX, VY , BVL , A AAA? 

50020FCRF  AT  SlH  ,6H THRUST ,E15.R,?X, 3HVX  ,t 15. 8,?X, 3HVY  ,fc 1 5. 8 ,2 X , 3HVL  , 
IE15.R,?X, 6HTHET A  ,E15.Bo 
PRINT  5003, AMASS, AX, AY, BL, AC 

50030FCRF  AT  SlH  ,6HMASS  , E 1 5. « , ?X , 3HAX  , E 15. B, ?X, 3HAY  ,E 1 5. 8 ,2 X, 3HAL  , 
ltl5.R,?X, 6HAG  , E 15.80 

PRINI  5004, AW, AT, AO, OEN, AM 

50040FCRMT  SlH  ,6HAW  , E15.H, ?X*  3UAT  , E 1 5. B, ?X , 3HA0  , E 1 5. 8 ,2 X , 3H0L N , 
IE15.B,?X, 6HAM  ,E15.8o 

PRIM  5005, SOS,  ANM,  At, D, ALT 

5C050FCRF  AT  SlH  ,6HS0S  , E 15. B, ?X, 3HNM  ,f 15.R,?X, 3HAE  ,E15.8,2X,3H  l<  , 
l  E 15 .8 « 2X , 6HALT  ,E15.Rn 
PRINT  5006, CRAG, AK l ,VW, OPW, GKAV 

50060F  CRN  AT  SlH  ,  6H0R  AG  ,  E  15.8,  ?X,  3HK  l  ,E  15.8,  2X,  3HVW  ,E  1 5.  8,2  X,  THOI’W , 
IE  15 .8 ,?X ,6HGRAV  ,El5.Rn 
IF  1SS5o701,701,910 
910  PUNCH  3900, T, BVL, HL, ALT ,X 
3900  FTRMT  S5F15.40 
701  I  0  I  £  CELT 

IF  SC  I linBO4,0O4, 796 

796  IF  SXXXX 087,797, 798 

797  XXXX  I  XXXX  £  1. 
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RMS  A  FRMAS 
THRUST  A  0. 

TAT-  CELT 
CELT  A  .1 
T  A  T  C  CELT 
BILL  A  0. 

TIRE  A  FTIME 
GC  TC  1903 

804  IF  «T  -  TIMED1903, 1903, 1900 
1900  TAT-  CELT 

CELT  A  TIME  -  T 
T  A  TIME 
BILL  A  BILL  Cl. 

1903  BIWT  A  BIWT  -t*RMAS«DELTn/32. 1740 
AMASS  A  BIWT 

IF  ISENSE  LIGHT  4u3901,3902 

3901  ALT  •  OVL  •  CELT 
GC  1L  3903 

3902  ALT  U  it  VY  *CELT 

3903  VL  A  BVL 

IF  il  -  T PO900 ,900,909 

900  IllEl  A  H  THETO 
AIP^A  A  0. 

GC  If  9720 

909  IF  ?T-tTP£AK 3 O 090 1, 90 1 , 872 
872  ALPT A  A  0. 

GC  1C  800 

901  ALPFA  A  AKA 

800  THE  1  A  A  AT  ANF3VY/VXH 
9720  IF  '.TIM!  -TrB7, 1400, 1400 
87  STOJ 
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FORTRAN  RUN 

HOP  16 

C  LAUNCH  PHASl  1 K AJECTORY  BY  W  B  WARRkN  FOR  RESEARCH  STAFF 

C  24  AUGUST  1962 

C 
C 

C  PROGRAM  NO.  649 

C  RMAS  »  PROPELLANT  MASS  RATE  OF  FLOW  -  LBS/SEC 

C  8IWT  A  TOTAL  INITIAL  WEIGHT  -POUNDS 

C  FWT  H  MISSILE  WEIGHT  AT  BURNOUT  -POUNDS 

C  VL  4  INITIAL  LONGITUDINAL  VELOCITY 

C  AL  »  MISSILE  LENGTH  -FEET 

C  A  H  VEHICLE  REFERENCE  AREA  -FEET 

C  TIME  t  TIME  OF  BURNOUT  -  SECONDS 

C  0  tt  DIAMETER  OF  MISSILE  -  FEET 

C  OTMETA  #  ANGLE  OF  THE  VELOCITY  VECTOR  WITH  RESPECT  TO  THE  HORIZONTA 
C  IL  -RACIANS 

C  AKl  »  THETA  AT  BURNOUT  -  RADIANS 

C  ALPHA  0  ANGLE  OF  ATTACK 

C  THRUST  0  POUNDS  FORCE 

C  AL T C  0  LAUNCH  ALTITUDE 

C  DR ACC  0  INITIAL  DRAG 

C  ALAM  0  LIMIT  ON  DRAG  COEFFICIENT  CURVE 

C  ALAM  0  LIMIT  ON  DRAG  COEFFICIENT  CURVE 

C  G  N  GRAVITATIONAL  FORCE  FOR  LAUNCH  ALTITUDE 

C  (HJNITS  Oi  TABLES  USED  -SPEED  OF  SOUND-F T/SEC ,  DENSI T Y-LB/F T»»J ,  GRA 
C  L  V  I  I  V  -FI/SFC*«2,  ALTITUDE  -FLLT 

C 
C 

flBHH  E  CRM  A  I  LL15.Hu 

RIAL  HHHR.SCLASS 

RE  AC  100, RMAS, BIWT, FWT, VL,Al 

RE AC  100, A, TIME, DEL ! , U, THETA 

KE  AC  100, AKl, AK 4, THRUST, ALTO, DR AGO 

REAL  100 , At  AM, G, AL  AMI, TP, »  T 1ME 

KLAL  100,  Al  A,  ERMAS 

REAL  100,SSl,SS?»SS3»SS4»SS5 

CALL  (LASS  LSCLASSu 

PRINT  1000, RMAS , RIWT , FWT ,Vl , AL 

PRIM  10n  | ,  A,  T  I  ME,  DLL  T  »  D,  I  HE  I A 

PR  INI  1002,  AKl  ,  AK4,  THRUST,  Al  TU,DI’.AG(I 

(.ALL  (IASS  hSf.L  ASSH 

1000  FORMI  X I ME , GMRMAS  , t 1 5. H, 2X , 6H I W I  ,E 15. 8,2X,6MF WT  ,E16.8,2X, 

I6MV1  ,Er».H,?X,ftH  L  ,115.811 

1001  THREAT  X  1 M  ,NH  A  ,  E 1 5. «, 2X, 6HT IME  ,£ 15. 8, 2X.6HOEL T  ,E15.8,2X, 

I6HDIAM  ,F15.H,2X,AHTHETA  ,El5.«u 

1002  TORFAl  l|H  ,6H  Kt  ,E15.R,2X»6H  K4  ,  E 15. 8,2X,6HTl<RUS!  ,E15.8,2X, 

1  Ml  Al  I  ('  ,E  1*1.8, 2X,6HDRAg0  ,115.8,//////////////////////////////// 

2  //////////////////////« 

SENSE  LIGHT  2 

XX XX  0  o. 

THI  10  0 T  HE  1  A 
AK i  0  TO. 


I 
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ALPHA  0  0. 

100  FCRMT  I5E15.8n 
T  0  0. 

CCCC  0  0. 

BILL  0  0. 

I  PAGE  0  0 
CCNTR  0  0. 

SENSE  LIGHT  4 
APASS  0  RIWT 
ALT  0  ALTO 

1400  CALL  TAP  IGRAV, ALT, ATEMK, APNMM, UEN, SOS.ATEMRo 
IF  *SSln7,7,6 

7  IF  ITHRUSTn  8,9,10 
9  SENSE  LIGHT  1 

GC  TC  7 3 

8  GC  TO  10 

10  THRU0THRUST 
GC  TO  ?3 

6  RE AC  INPUT  TAPE  6, 2100, T l, THRU1 
2100  FORPAT  I2E15.8n 

IF  *T-Tla7200, 7201, 7202 

7201  THRU  0  THRU1 
THRU2  0  THRU1 
GC  TO  23 

7200  THRU20THRU1 
T2  0  Tl 
GC  TO  6 

7202  TRAT  0»T-T2o/tTl-T2n 

THRU  0  THRU 2  £  IITHRU2-THRU la*  TR ATd 
BACKSPACE  6 
GC  TO  23 
23  AM VL/SCS 

I F tAPa8  7,76,11 
76  ORAC0O. 

GO  TC  151 

11  CALL  ORAGCO  tORAG, AMACH, AMa 
151  IF  ISS2o27, 27, 26 

27  OPV.  00. 

GO  TO  2222 

26  CALL  WIND  (ALT , VMa 

31  OPw0IOEN»VW**2o/2. 

2222  AO»3«l./2.a»DEN*<VL**2a*ORAG«Aa/AMASS 
IF  IT  -  TPal53, 153,915 
915  IF  ISS3o777, 777, 32 
777  IF  !ALPhAd152, 153, 153 

32  REAC  INPUT  TAPE  7, 2300, T IM 1 , THE T  1 
2300  FORPAT  12E15.8a 

IF4T-TIP1034, 35,36 

35  THri2  0THET 1 
GC  TO  152 

36  TIP2  IT  INI 
THET20THET1 
GO  TC  32 
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34  TRATAtT-T I M 2a/ XT IMl-T IM2u 

THE  T2  ATH£T2-XTRAT*XTHET?-THETIcd 

152  REAC  INPUT  TAPE  9, 2 300. T !M 3, ALPH 1 
IF  XT-TIM3u37.38.39 

38  ALPHA  A ALPH1 
GO  TO  153 

39  TIM4ATIM3 
ALPH2  IT  ALPH  1 
GO  TO  152 

37  T3RAT  AXT-T  IM40/XT  IM3-T  lM4n 

ALPHA  A  ALPH2  -XT 3R AT*X ALPH2-ALPHlnn 

153  AE  A  AL  *0*S INFXTHET  Aa 
AMAXCPW*AEa/ AMASS 

IF  XSENSE  LIGHT  la  700.701 
700  IF  XSS4a  702.702.703 

703  REAC  INPUT  TAPE  8. 2200, FT IME 1 , ALA l 
2200  FORMAT  X2E15.8o 

IF  tT-FTlMEa704.705.706 
706  FTIME2  AFT  IME 1 
ALA?  A ALAI 
GO  TO  703 
705  ALA  A ALA l 
ALA2A  ALAI 
FT  I ME2AF1 IME1 
T  A  F  T I  ME  1 
r.C  IC  70? 

704  TIMtR  A  TT-FUME20/XFTIME1  -  FT IMC2o 
ALA  AALA?  -tT IMER*XFT IME2  -FT IME loo 
BACKSPACE  8 

70?  If  tCCr.r.n769,769,770 

769  KEAL  100, THRU, TTTT 
CCCC  A  CCCC  £  l. 

GO  TC  911 

770  IF  tl T TT n50 1,501,911 

501  THRU*  XAMASS/CnSFXALPHAau*XAL ACADCAG*  SINFXTHE TAa-AM*COSFXTHETAau 
911  AT  A THRU/ AMASS 

AX  A  AT  *COSF£THET  A  -  ALPHAn  -  X AD*COSFX  THE T Ann£  AM 
IF  TAX  -  . lo925.925.926 

925  AX  If  0. 

926  AY  A  AT  *S INFXTHET  A  -  ALPHAn  -  XAD*SINF  XTHE TAaa-GRA V 
RL  A  AX •  COSFXT  MET  An  C  AY*  SINUTHETAo 

IF  tSENSl  LIGHT  4n51,52 
51  VX  AO. 

VY  A  0. 

BVLA  HL 
X  A  0. 

Y  A  0. 

CVL  A  0. 

SET  SI  LIGHI  4 
GC  1C.  54 

52  VX  A VX  L  XAX  •  CELT  o 
VY  A  VY  f.  XAY  •  OF.LTn 
902  BVI.  A  SCRT FXVX **2  £  VY**2o 
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X  A  X  E  VX*C£LT 
Y  I  Y  £  VV  •  CELT 
CVL  A  SCRTF*X«*2  €  Y«*2o 
54  AG  A  GRAV 

296  CGNTR  ACCNTR  £  I. 

IF  *CCNTR-l.o297, 298, 299 

297  PAUSE 

299  IF  tCONTR  -7.Q399, 399, 365 
365  CGNTR  A  0. 

PRINT  9462 
9462  FORNAT  XlHKn 

CALL  CLASS  tSCLASSo 
GO  TC  296 

298  (PAGE  A  1  PAGE  C  I 
PRINT  5000, IPAGE 

5000  FGRMT  X  1H1 ,  11  OX,  5HPAGE  ,  I  ln 
CALL  CLASS  XSCLASSn 
399  PRINT  4999 
4999  FORMAT  X1H  ,//D 
5701  AAA A |  A  ALPHA/. 01745 
AA AA2  A  THETA/. 01745 
PRINT  5001, T,X,Y, CVL, AAAA1 

50010FDRMT  XlH  ,6HT  1ME  ,E15.R,2X,3H  X  ,E15.8,2X,3H  Y  ,E15 
1E15.8.2X.6HALPHA  ,E15.8o 
ANN  1  X/6080. 

ALT  A  Y  C  ALTO 

PRINT  5002, THRU, VX, VY, BVL , AAAA2 
50020FGRKAT  XlH  ,6H THRUST , E 15.8, ?X, 3HVX  ,6 t 5. 8, 2X, 3HVY  ,E13 
1E15.8, 2X.6HTHETA  ,E15.8o 
PRINT  5003, AMASS, AX, AY, BL, AG 

50030FGRMT  XlH  ,  6HMASS  , E15.8,2X,  3HAX  , E  1 5.  0,  2 X ,  3HAY  ,E15 
IE15.8.2X.6HAG  ,E15.8o 

PRINT  5004, AW, AT, AO, OEN, AM 

50040FQRN  AT  XlH  ,6HAW  ,615.8, 2X, 3HA T  , 11 5 . 8 , ? X , 3HAD  ,E15 
1E15.8.2X, 6H AM  ,E15.8u 

PRINT  5005, SOS, ANM,Afc,D, ALT 

50050FORMT  till  ,6HS0S  ,  E  15.8,  2X,  3HNM  , l:  1  5.  8,  2X,  3HAE  ,E15 
lfcl5.8,2X,6HALT  , E  15.80 
PRINT  5006, CRAG, AK1,VW,DPW, GRAV 
50060FCRI AT  XlH  , 6H0R AG  , 6 1 5 . 8, 2X , 3HK 1  ,F 15.8, 2X, 3HVW  ,E15 
IE15.8.2X , 6HGRAV  ,E15.8n 
IF  *SS5u70l , 701,910 
910  PUNCH  3900, T, BVL,  BL  , ALT, X 
3900  FORM!  X5F15.4U 
701  T  A  T  £  UFLT 

If  .111110804,804,796 

796  IF  XXXXXU87 , 797* 798 

797  XXXN  A  XXXX  £  1. 

RPA*.  A  FKMAS 
THRUST  A  0. 

TAT-  OUT 

Cl  L 1  r  .1 
I  A  T  t  CELT 


8.2X.3H  L  , 

8.2X.3HVL  , 

8,2  X , LHAL  , 

8 , 2  X , 3H0E  N , 

8 ,2 X , 3H  U  , 

8 , 2  X , 3HDPW » 
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BILL  It  0. 

TlHf  H  FT! HE 
GC  TO  1903 

804  IF  SIT  -  T  IM£nl903, 1903,  1900 
1900  T  H  T  -  CELT 

DtLI  H  T IHE  -  T 
T  H  1 IME 
BILL  H  6ILL  Cl. 

1903  BIWT  H  BIWT  -URMAS*DEL To/ 32 . I 74Q 
AP ASS  U  0 IWT 

IF  S1SENS F.  LIGHT  4a3901,3902 

3901  ALT  tt  BVL  *  CELT 
GC  TC  3903 

3902  ALT  It  YC  V Y  -CELT 

3903  VL  It  BVL 

IF  *T  -  T  Pn900 , 900, 909 

900  THETA  U  THETO 
ALPHA  It  0. 

GO  TC  9720 

909  IF  tT—?TPEAK3uo901, 901,872 
872  ALP!  A  H  0. 

GC  TC  800 

901  ALPHA  It  AK4 

800  THE  1  A  A  AT  ANFIVY/VXo 
9720  If  *UFE  -Tnfi7, 1400,  1400 
87  STOP 
798  STOP 
ENC 

BOP  CLASS  16 

PRUGRAPPEC  BY  W  8  WARREN 
CLASSIFICATION  OF  WORKING  PAPERS 
SUBROUTINE  CLASS  *SCLASSo 
IF  ISCLASSO4000, 6000,6001 

6001  IF  1SCLASS  -  l. 04000,6002, 6003 

6003  IF  UCLASS  -  2 .04000,  6004, 4C00 
6000  °R  I  AT  6010 

6010  FCRPAT  *1H  *  35X , 12HUNCL  ASS  I F I  EDO 
PRIAT  6011 

6011  FORPAT  *1H  , 34X , 14HW0RK ING  PAPERSo 
GC  IC  6013 

6002  PRIAT  6012 

6012  FORPAT  *1H  , 32X ,  12HC0NF  IDEN  T  IAL o 
PRIAT  6011 

GC  TC  6013 

6004  PRIAT  6014 

6014  FCRMT  tlH  ,  35X,  6HSkCRETo 
PRIAT  6011 
GC  TC  6013 
4000  STOP 

6013  RETURN 
END 
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BOP  TAP  16 

PROGRAMED  BY  W  B  WARREN 
TRAJECTORY  ATMOSPHERIC  PROPLRTItS 

SUBROUTINE  TAP  TGR AV t AL T, AT EMK , APNMM, DEN , SO S. A TEMRn 
IF  TSENSE  LIGHT  6D1  l, 13 

11  OIMEKSICN  APT9.6D 
SENSE  LIGHT  A 

00  12  I  A  1,9 

12  RE  AC  105, APfI,in,AP*I,2o,APTI,3c,AP<I,6o 
105  FORMAT  T6E15.8D 

RE  AC  100, UMTS 
100  FORMATTE 15.8a 

13  IF  IUMTSD500,  1000,  500 
1000  ALT  A  ALT/3. 28083989 

VL  A  VL/ 3.28083989 

500  GR AV  A9.80665-T6371239.9/T6371239.9  E  ALTOO--2 
H  A  ALT-6371239. 9/T6371239. 9  C  AL  Tn 

DO  202  I  A  2,9 

IF  t H  —  APT  I , loa203,202,?02 

202  CONTINUE 
I  A  10 

203  J  A  !-l 

ATEMK  AAP«J,2ut*H-AP«J,  lOD-APTJ,  3D 
IF  TAPTJ, 3no206, 205,206 

206  APNMMA APTJ, 6n-TAPTJ,2n/ AT GMKU--T. 036 1 6679/APt J,3nu 
GO  TC  206 

205  APNPMAAPTJ, 6n/ ETPFTT.Q36 166  796D-XH- APT J , laa/APTJ,2na 

206  DEN  AT*  .00368 38 39 AD -APNMM /A  TEMKd 
SCS  A20.0A6  333-AT FMK- • . 5 

I FTUM  TSo50 1 , 502*501 
502  GR AV  A  GRAV -3. 28083989 
HAH-  3.2808389 
ALT  H  ALT  •  3.28083989 

CEN  A  OtN  •  XX2.205/X3.280839B9--3DD/3?. 176a 

SCS  J  SCS  •  3.28283989 

501  ATLMH  A  TAT  EMK— 273. 16D-T9. / 5.u  L  691.69 

IF  3 AL T  -  60J000.D506.506, 505 

505  DEN  A  0. 

SCS  AO. 

506  HETUKu 

em: 

8UP  WINO  16 

PROGRAMMED  BY  W  B  WAKREN 
WlNt  I'KCFILE 

SUHHOtlT  INE  WIND  <ALT,VWO 
IF  TSi-NSC  LIGHT  6nU,13 

11  OIPC.V.ICN  AWTI  3,  3d 
SENSE  LIGHT  6 

DC  12  111,13 

12  REAL  105,  AWtl,  la,  AWTI,  2a,  AWT  I,  3c 
105  FCKPAT  X 3E 15.8a 

13  DC  202  112,13 
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IF  HALT  -  AWX1, lna203,202, 202 

202  CONTINUE 
I  H  13 

203  J  »  l-l 

VU  #  AWHJ.2Q  £  HALT  -  AWX J ,  lnu*AWXJ , 3n 
SOI  RETURN 
END 


BOP  OR AGCO  16 

C  PROGRAMMED  BY  W  B  WARREN 

C  ORAG  COEFFICIENT 

SUBROUTINE  CRAGCO  XCRAG. AMACH, AMo 
IF  XSENSE  LIGHT  2011,13 

11  DIMENSION  ACUGX9,4o 
DO  12  K81.9 

12  READ  105 « ADRAGXK,  la, ADR AGXK  ,  2n*  ADRAGXK ,  3a « ADRAGXK  ,4a 
105  FORMAT  X4E15.8a 

13  00  20?  K«2.9 

IF  HAM-ADRAGXK, loo203,202,202 

202  CONTINUE 
K  »  9 

203  J  A  K-l 

IF  XAORAGXJ ,4aa204,  205,204 
205  ORAG  8  AM  *  ACRAGXJ.2D  £  ADRAGXJ.3D 
GO  TO  350 

204  IF  HAC.RAGXJt  In  - 1 .  Id300,  300.  301 

300  ORAG  H  XXXAM-ADRAGXJ, 2ao**2D/A0RAGXJ,4nn  £  A0RAGXJ,3a 

GO  TO  350 

301  UHAG  U  A0RAGXJ.3D  -  XSORT FX AORAGXJ , 4o*XAM-A0R AGtJ ,2ooaa 
350  AMACH  A  AM 

500  RETURN 
ENC 


